Combinatorial nuclear level-density model 
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A microscopic nuclear level-density model is presented. The model is a completely combinatorial 
(micro-canonical) model based on the folded- Yukawa single-particle potential and includes explicit 
treatment of pairing, rotational and vibrational states. The microscopic character of all states 
enables extraction of level distribution functions with respect to pairing gaps, parity and angular 
momentum. The results of the model are compared to available experimental data: neutron separa- 
tion energy level spacings, data on total level-density functions from the Oslo method, cumulative 
level densities from low-lying discrete states, and data on parity ratios. 



I. INTRODUCTION 

Nuclear many-body level-density models are key ingre- 
dients in nuclear reaction theories, where they, for exam- 
ple, govern the rates and decay patterns of astrophysi- 
cal processes and nuclear fission. In statistical methods, 
for example the Hauser-Feshbach formalism [l| for de- 
scribing nuclear reactions, a knowledge of the level den- 
sity is crucial @, d, 0|- How to calculate the nuclear 
level density (NLD) has been a long-standing challenge 
[E B 0i B S EH • Recently it has been subject to renewed 
interest, theoretically as well as experimentally. 

The simplest type of model is the Fermi-gas model, 
which is based on the partition-function method. It pro- 
vides simple analytical formulas for the NLD [ll[ . Several 
phenomenological extensions have been proposed in order 
to reproduce experimental data. By adjusting free model 
parameters to data these models give unprecedented ac- 
curacy in the region of the parameter fit 0, [13] . Also 
semi-classical methods have been used to obtain expres- 
sions for the level density [H|. However, the Fermi-gas 
models are unreliable outside these regions, eg. when 
they are extrapolated to higher excitation energies or 
to nucleon numbers far from stability. Ideally nuclear 
structure should be included in NLD models. Several 
combinatorial models based on nuclear mean-field the- 
ory have been proposed, see eg. Refs. OH!]- Beyond 
mean-field methods have also been used to model the 
NLD, eg. the Shell-Model Monte-Carlo method [H [T^ 
and the interacting shell model [l8| . These latter models 
take into account effective nucleon-nucleon interactions, 
but at the same time suffer from limitations due to the 
limited size of the Hilbert space, and hence are presently 
unable to provide global predictions for level spacings at 
the neutron-separation energy. 

Experimentally the NLD has been subject to renewed 
interest in the last decade partly due to the development 
of the Oslo method, which has provided new experimen- 
tal data [l| . The Oslo method provides the level density 
over extended regions of excitation energy as opposed 
to the neutron-separation level-spacings data which only 
provide one data point at relatively high excitation en- 
ergy. Also recent measurements of separate level densi- 



ties of 2 + and 2 states [20| challenge theory to reproduce 
these observed parity ratios. 

Few nuclear-structure models have been used to si- 
multaneously globally describe nuclear masses, fission 
barriers, ground-state spins and decay rates. One such 
model is the microscopic-macroscopic FRLDM model 
which has previously been used to model these observ- 
ables [H, [U [12, HH , and here serve as the starting point 
for calculating the nuclear level density. In this paper 
a combinatorial (micro-canonical) nuclear level-density 
model based on the folded- Yukawa single-particle model 
is presented. The model is fully microscopic with pairing 
correlations, vibrations, and rotational excitations cal- 
culated for each many-particle-many-hole excited state. 
Presently, the lowest ten million or so states can be ac- 
counted for. This implies that the excitation energy re- 
gion from the ground state to well above the neutron res- 
onance region is included. No additional parameters are 
introduced in the model, and no refitting of parameters 
of the FRLDM is performed. 

In Sec. |TT] the combinatorial folded- Yukawa (CFY) 
level-density model is described. The model allows ex- 
plicit tracking of quantum numbers, and distributions 
of pairing gaps, parity and angular momentum are dis- 
cussed in Sec. Mil Results from the CFY model are com- 
pared to experimental data in Sec. IIV1 and in Sec. fVlthe 
CFY model is compared to other theoretical NLD mod- 
els. Finally, a short summary is given in Sec. I VII 



II. 



THE COMBINATORIAL MODEL OF NLD 
(CFY) 



The NLD is calculated by means of a combinatorial 
counting of excited many-particle-many- hole states as de- 
scribed in Sec. Ill Al In Sec. Ill B I we present how pairing 
is taken into account for excited states by explicitly solv- 
ing the BCS equations for all individual configurations. 
Rotations are taken into account combinatorially with 
a pairing-dependent moment of inertia (Sec. Ill Cp . The 
vibrational contribution to the NLD is investigated by 
including microscopically described phonons using the 
Quasi-particle Tamm-Dancoff Approximation (QTDA), 
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see Sec. lII Dl In Sec. Ill El it is described how we briefly ac- 
count for a general residual interaction causing a smear- 
ing of the level-density distribution at higher excitation 
energies. 

All produced nuclear levels, calculated as described be- 
low, are sorted into a binned level density, where the typ- 
ical bin size is in the range AE = 30-50 keV. The level 
density is calculated by counting the number of levels in 
the energy bin, 

p(E b ,I,n) = -^= { " ' ^SiE-Eiil^^dE, (1) 

where Eb is the energy center of bin b and Ei (J, n) 
denotes the calculated state with energy Ei (given by 
Eq. (|15p ). angular momentum / and parity tt. The total 
level density at a given excitation energy E is 

p tot (E)=J2p(EJ^). (2) 

J,7T 

A. Combinatorial intrinsic level density 

A high-quality combinatorial level-density model re- 
quires a realistic description of the single-particle ener- 
gies. We shall here utilize a well-tested mean-field model, 
namely the microscopic-macroscopic finite range liquid 
drop model (FRLDM) [2l|, that provides a good global 
description of several nuclear-structure properties such as 
masses, fission barriers and beta-decay properties. The 
good agreement between ground-state spins calculated 
in the FRLDM model and experiment implies that the 
single-particle spectrum close to the Fermi surface is well 
described [22| . 

The single-particle energies are thus obtained by solv- 
ing the one-body Schrodinger equation, 

{T + V FY {e))\v)=e v \v), (3) 

for protons as well as for neutrons, where Vfy is the 
Folded- Yukawa single-particle potential. The parameters 
of the potential, as well as the deformation parameters, 
e = {e^i £3, £4, ■•■), are taken from an extensive calcula- 
tion of nuclear masses [2l| . All parameters of the single- 
particle model are thus fixed from other studies. One 
aim of the present study is to see how well highly excited 
states (in terms of level-densities) are described for all 
nuclei heavier than 16 O. 

The many-body ground state (the many-body particle 
vacuum) |0), is obtained by filling the lowest N (Z) states 
(including time-reversed states), 

N (Z) 

10} = £ a+|->, (4) 
v=i 

where |— ) is the one-body vacuum. Intrinsic excitations 
are obtained by many-particle-many-hole excitations on 



the many-body particle vacuum, 

n 

K)=n< a ^i°)' ( 5 ) 

Q = l 

where v and v' span all single-particle states (as well as 
time-reversed states) in the potential. We include all n- 
particle-n-hole states with n < 10 for protons as well as 
for neutrons. Corresponding energies are given by 

n 

E i = E + J2( e »o,- e K S )> ( 6 ) 

Q = l 

where Eq is the ground-state energy. (When pairing 
is included this expression is modified as described in 
Sec. Ill Bl ) Nuclei are assumed to have constant defor- 
mation for all excitation energies. For spherical and 
well-deformed nuclei this assumption is expected to be 
a good approximation for excitation energies here con- 
sidered (mainly below the neutron separation energy). 
For transitional nuclei this approximation might, how- 
ever, introduce some inaccuracies in the calculations. 

In addition to the energy and the number of unpaired 
nucleons (seniority), we also register the total parity and 
-fsf-quantum number for each excited state, where K as 
usual is the angular momentum projection on an intrinsic 
symmetry axis. The identification of a if-quantum num- 
ber can be done for all nuclei by assuming an infinites- 
imal deformation for spherical nuclei. The i^-quantum 
number of the individual unpaired nucleons are allowed 
to couple in all different ways to build up the total K- 
quantum number. 



B. Pairing 

The many-body wave-function of the excited states is 
approximated by the BCS wave-function with excited 
quasi-particles 

V) = n (-k» + ev4,,4„)x 

f"£T 2 

x n a i n V* + v » a t a t) i°> . (7) 

I/'GTl V&To 

where T2,t\ and To denote the spaces of double, single 
and zero quasi-particle excitations, respectively, used to 
build the many-body state i of Eq. [5] U v and V v are the 
standard BCS vacancy and occupation factors and |0) is 
the particle vacuum, see e.g. Ref. (24|. For the excited 
pairs in the group T2 the effect (compared to zero-quasi- 
particle states) is simply 

U v -> -V v , V u -» Uv (8) 

The pairing gap A and the Fermi energy A are obtained 
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by solving the BCS-equations 



Rotations 



A = G 



n = 2J2 v » + J2 1 + 2 J2 u *> 



(9) 



(10) 



VETO 



for each state. 

The pairing strength in the BCS model is governed 
by a single, primary BCS model parameter, namely the 
parameter of the effective-interaction pairing gap (r in 
Eq.(47) in [25[ or the equivalent c p in Eq. (49) and 
c n in Eq. (50) in [26]). In an adjustment of macro- 
scopic, pairing, and other model parameters to optimize 
a nuclear mass model r — 4.8 MeV was obtained for 
the standard BCS pairing model [2l|. Because the BCS 
model we use here, in contrast to Ref. [2l|, [2||, includes 
blocking and particle-number-projection it is necessary 
to use a pairing strength optimized for this method. In 
Ref. [26[ c p = c n = 4.95 MeV were obtained in an ad- 
justment of calculated pairing gaps to odd-even mass 
differences. However, it has been pointed out [U [25| 
that odd-even mass differences are subject to numerous 
non-smooth contributions, for example from deformation 
changes and from irregularities in the microscopic level 
structure, not just odd-even staggering due to pairing. It 
is therefore better to determine the pairing strength from 
a full nuclear mass calculation that includes the particu- 
lar pairing model under consideration and the associated 
adjustment of all model parameters [21(. In the folded- 
Yukawa macroscopic-microscopic mass model this yields 
the value 4.5 MeV for the optimized strength parameter 
of the particle- number-projected BCS model 26]. The 
pairing strength G, which is used to calculate pairing 
gaps for all excited states in the level-density calculation, 
is determined from this average pairing gap for each nu- 
clear system; for details see Ref. 26]. 

It follows that the excitation energy of the intrinsic 
many-body configuration for one nucleon type is 

E mh , t = 2 e » V ? + X e - + 2 X e » U l<<- 

VETQ V'ETl !/"ET2 
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where t denotes protons or neutrons and E® is the proton 
or neutron part of the ground-state energy. The pairing 
gap and Fermi level are calculated by solving the BCS 
equations Eqs. ([9]) and ([TO]) . The total intrinsic excitation 
energy of a many-body state i is 



E 



mb 



E, 



mb,p 



E, 



mb.i 



(12) 



The combinatorial approach to calculate the level den- 
sity involves calculating all possible many-body states, 
B^ b , which thus means that the BCS equations are 
solved about 10 7 times for each nucleus. 



A general feature of deformed nuclei is the existence of 
rotational bands built on the ground state as well as on 
excited states. The question is how to incorporate these 
states in the NLD. In principle, rotational states may be 
microscopically treated by solving the cranking Hamil- 
tonian, h u = h — uJrotjx, where uj rot is the rotational 
frequency and j x is the angular momentum operator for 
rotations around the x-axis. Matrix elements of the j x - 
operator correspond to energy excitations of the order 
AE = E2^uj « 62 ■ 41j4~ 1 / 3 MeV for a nucleus with mass 
number A and quadrupole deformation Ei [27| . This en- 
ergy can be compared to a typical energy of lowest ro- 
tational state, E 2+ = |^2(2 + 1) w 90^ 5 / 3 MeV, and 
we see that AE >> E 2 + for all considered nuclei. Since 
the collective state and the corresponding matrix element 
are well separated in energy, the risk of double-counting 
states by adding a rotational band on a microscopically 
calculated band- head is quite small, cf Ref. This is 
true as long as the excitation energies are smaller than 
energies corresponding to the temperature, T — AE, i.e. 
for excitation energies smaller than about 10A 1 / 3 MeV 
(for £2 — 0.25), giving about 30 MeV for mass A=25 
and 55 MeV for A=160. For higher excitation energies a 
saturation of the rotational enhancement should set in. 
In this study we concentrate on lower energies, say up to 
about 10 MeV. It is therefore a reasonable approximation 
to account for rotational enhancement by simply adding 
a rotational band on each band head; double-counting of 
states would not occur. 

Rotational states are consequently taken into account 
by adding a rotational band on top of each intrinsic band- 
head for deformed nuclei (here defined as nuclei with cal- 
culated quadrupole deformation | £2 1 > 0.05). The rota- 
tional energy of the different angular momentum states, 
J, in the rotational band is given by 



I(I + 1)-K? 



(13) 



2Jl(e a ,A*,Ai)' 

where Ki is the angular momentum projection on the 
symmetry axis of the intrinsic state i upon which the 
rotational band is built. The quadrupole deformation is 
denoted by e 2 while A p and A„ denote the proton and 
neutron pairing gaps of the intrinsic state i. The moment 
of inertia around an axis perpendicular to the symmetry 
axis i7j_( £ 2> A* A*) is approximated by the rigid-body 
moment of inertia with deformation 62, modified by the 
calculated pairing gaps for the considered state, as given 
in Ref. [281 ]. Given the angular momentum projection 
K and parity 7r of the band-head the rotational band 
includes the following levels 



K W ,(K + l) 7r ,(K + 2)V- if if 7^0, 
I n = ( 0+,2+,4+,... if K = 0^ 

l-,3-,5-,... if -fiT = 0" 



(14) 



The Coriolis anti-pairing effect is neglected and no vir- 
tual crossings of rotational bands are taken into account. 
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FIG. 1: (Color online) Rotational enhancement, K to t, for the 
nucleus 162 Dy as a function of excitation energy. The insert 
shows the enhancement compared to the simple enhancement 
model of Ref. [29( ] (red dashed line). 



Thus, the pairing gap and moment of inertia are assumed 
to be unchanged from the band-head pairing gap for all 
states in the rotational band. This approximation is rea- 
sonable since mainly low-spin states play a role in the 
present study. 



by 



The total energy of the calculated level is thus given 



Ei=E. 



nib 



(15) 



where E^ h is the energy, Eq. (fT2")l . of the intrinsic many- 
body configuration i, and El ot {I,Ki) is the rotation en- 
ergy, Eq. I|13p. of the level with angular momentum I 
built from the intrinsic many-body configuration with 
angular momentum projection Ki. These energies are 
used to calculate the level density according to Eq. (fTJ) . 

Fig-Hshows the rotational enhancement, K Iot , for the 
well-deformed nucleus (e 2 =0.26) 162 Dy (cf. Fig. [T2J cal- 
culated as the ratio of the level density when rotations 
are included or excluded. For low excitation energies 
there are large fluctuations which are artifacts of the low 
level density combined with the smoothing procedure of 
Sec. Ill El For higher excitation energies (> 3 MeV) the 
rotational enhancement is a slowly increasing function, 
of the order of a factor 5 at the neutron separation en- 
ergy. This prediction is compared to the SU(3) model 
of Ref. [1^1, which is shown in the insert of Fig. [1] The 
SU(3) model gives almost an order of magnitude larger 
enhancement for excitation energies in the region of the 
neutron separation energy. We note that a combinatorial 
level-density model based on the Nilsson potential gives 
a rotational enhancement for 162 Dy similar to our results 

M. 



D. Vibrations 

Many nuclei exhibit low-lying states of vibrational 
character, which are usually of quadrupole or octupole 
type. Such low-lying vibrational states appear at exci- 
tation energies not much lower than the 2qp excitations 
which in a coherent way build up the collective state. We 
therefore believe that the vibrational enhancement of the 
level density, contrary to the rotational enhancement dis- 
cussed above, must be described microscopically. 

In order to describe vibrational states the Quasi- 
particle-Tamm-Dancoff- Approximation (QTDA) is used. 
According to the Brink- Axel hypothesis [3l|, [32| phonons 
are built on every intrinsic many-body configuration 
^mb- T ne QTDA equation is solved for each state i in 
order to get all possible phonon excitation energies and 
wave- functions. This means solving the QTDA equations 
millions of times for every single nucleus. 

The residual interaction is approximated by the dou- 
ble stretched (isoscalar) Quadrupole-Quadrupole inter- 
action. This interaction is well defined in the case of a 
harmonic oscillator potential. In the case of a finite-depth 
potential as the folded- Yukawa potential the interaction 
should take into account additional finite size effects, for 
example as is done in Ref. (33[. In the present work the 
finite-depth effects are ignored and the double stretched 
approach is used as defined in Refs. [34l l35l] . 

The QTDA secular equation can be written [3(| 



1 

X2K 



(16) 



where the effect of Eq. ([8]) has not been explicitly written 
out. Qik is the double stretched quadrupole operator, 
where the components K = and K = 2 are considered. 
The set of roots {(Tiu;)*-} of this equation is the excitation 
energies of the vibrational phonons j on top of the intrin- 
sic state i, whereas the poles E^ 9 ' 1 = \J (e^ — A) 2 + A 2 
are the unperturbed two-quasiparticle excitations on the 
many-body configuration E l mh with pairing gap A, as 
calculated by Eqs. (|9|) and (fTf]| . and e M are the single- 
particle energies. 

The self-consistent coupling strength is given by UK 



X2K 



8tt 



Muj 2 



Ai 



92K 



A(Q m )' 



(17) 



where 520 = 1 and c/22 = — 1. The expectation- values 
(f 2 ^ and (Q20) are calculated in double stretched co- 
ordinates [3J|. To test that this gives overall reason- 
able result, we have verified for a large number of nuclei 
that for the K = and the K = 2 components of the 
iso-scalar Giant Quadrupole Resonances the calculations 
agree well with the systematics of the Giant Resonance 
energies ftw G QR = MeV [ai]. 

The phonons are never repeated and double counting 
is explicitly avoided by the following procedure. The 
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phonon wave functions are given by 



4a„ 



(18) 



where X^ jV are the wave-function components of all ex- 
cited quasi-particle states aj^a^ on top of the configura- 
tion E^ h . The level density is increased by one state at 
the energy of the phonon (Tkj)J-, and decreased by the 
amount given by the wave-function component X 2 v at 
the energy of the corresponding pole (E^ p ' 1 +E^ P,Z ). The 
change in level density due to the phonon j on top of in- 
trinsic state i is thus 



5p{E) = 5(E - E. 



mb 



(Mj)- 



E-K 



mb 



ET 1 )), (19) 



where E is the excitation energy relative to the ground- 
state. This procedure increases the level density by re- 
ducing the strength of the pure quasi-particle excitations 
and adding strength by means of the phonons (which are 
pushed down in energy due to the QQ-interaction). 

Most phonons have little collectivity in the QTDA ap- 
proximation, as most of the phonon excitations lie very 
close to a pure quasi-particle excitation, and hence the 
wave-function is completely dominated by that single 
component. Even if there are a few very collective low- 
lying phonon states the vibrational enhancement of the 
level density becomes small, since most excited states 
provide quite non-collective phonon excitations. Due to 
this non-collectivity of most phonons there is no way that 
the phonons, in general, can be repeated to form two- or 
even thrcc-phonon excitations. 

The vibrational enhancement factor in this method is 
in general quite small, of the order of a few percent. This 
is substantially lower than predictions from other models. 
For example, the attenuated phonon method [13. ILH |3?]] 
and the Boson partition function method [38[ both give 
up to an order of magnitude enhancement at the neu- 
tron separation energy. Fig. [2] shows the vibrational en- 
hancement as a function of excitation energy for 162 Dy. 
The effect is very small, close to 1 % at 7 MeV excita- 
tion energy. For the same nucleus the attenuated phonon 
method gives an enhancement factor of about 3 as shown 
in the insert of Fig. O 

The level-density enhancement due to quadrupole vi- 
brations is thus found to be very small. Higher multipolc 
vibrations (as octupole vibrations) are expected to con- 
tribute with enhancements of the same order of magni- 
tude or less and are neglected in the CFY model. 
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FIG. 2: (Color online) Vibrational enhancement, if v ib, us- 
ing the QTDA method for the nucleus 162 Dy as a function 
of excitation energy. The insert shows the enhancement com- 
pared to the enhancement of the attenuated phonon model of 
Ref. (red dashed line). 



from the combinatorics. Smearing effects from the resid- 
ual interaction can approximately be taken into account 
by assuming a spreading width of all excited states. The 
spreading width is implemented in terms of a Gaussian 
envelope with width a, i.e. the delta functions in Eq. [T] 
are replaced by Gaussians. Estimates of the spreading 
width FWHM gives [H 



/ A \~ 1/2 
r = 0.039 ( — J E 3/2 MeV, 



(20) 



where the FWHM is related to the Gaussian spreading 
width (7 



2V2 In 2 ' 

Assuming that all excited many-body states in an en- 
ergy bin AE are uniformly distributed, the level density 
becomes 



p{E b ) = Y J P{Ea)\ 



erf 
erf 



E a + AE/2 - E b 

E a - AE/2 - E b 
V2~a 



, (21) 



for bin-point b. The method implies a smearing out of 
level-density properties over a range T, which is smoothly 
increasing with excitation energy. As a consequence fluc- 
tuations of energies and wave- functions in the range T 
follow GOE statistics of Random Matrix theory, that is 
often denoted as quantum chaos in the nucleus, see e.g 
Refs. IB 



E. Many-body damping width 

In the mean-field approach all excited many-body 
states are treated as non-interacting. A residual two- 
body interaction will mix the many-body states obtained 



III. LEVEL DISTRIBUTIONS 

In the present combinatorial approach it is possible 
to extract different distributions exhibiting details of the 
level-density function. We present here microscopically 
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calculated distributions for pairing gaps (Sec. lIII A")) , par- 
ity (Sec. IIIIBj) . and angular momentum (Sec. IIIIC]) . 

A. Pairing-gap distribution 

The BCS equations, Eqs. ([9]) and (fTT)|) . are solved for 
all individual many-body configurations and hence pair- 
ing gaps for all states are obtained. In Fig. [3] the dis- 
tribution of proton pairing gaps for 162 Dy is shown for 
a number of excitation energies. For the lowest excita- 
tion energies only the ground-state and the states in the 
ground-state rotational band exist. The pairing gap of 
the levels in the rotational band is fixed to be the same 
as the pairing gap of the band-head, see Sec. Ill CI As the 
excitation energy increases levels with reduced pairing 
gaps appear. However, no transition to a completely un- 
paired system is observed, and levels with non-collapsed 
gaps (A > 0) survive to the highest considered excitation 
energies. At the highest considered excitation energy for 
162 Dy, £ oxc = 8.4 MeV, 36% of the levels still have a non- 
zero proton pairing gap, see Fig.[3j The mean value of the 
proton pairing gap (A) at different excitation energies is 
shown in Fig.Hfor 162 Dy. Between 2 MeV and 3.5 MeV 
excitation energy there is a rapid decrease in the mean 
pairing gap. At higher excitation energies the decrease 
is slower. At 8.4 MeV excitation energy the mean value 
is (A) = 0.2 MeV. The reason why pairing may survive 
and be of substantial size in states at these high ener- 
gies can be explained as follows: Several highly excited 
states are built by exciting particles and holes far from 
the Fermi energy, where blocking in the BCS equations 
has no effect on the pairing gap. 

The non-collapsed pairing gaps influence the moment 
of inertia and keep it reduced as compared to the rigid- 
body value. This implies that even at excitation energies 
in the region of the neutron separation energy the mo- 
ment of inertia is on average smaller than the rigid-body 
value. A similar observation is made in Ref. 11511. 



B. Parity distribution 

In Fermi-gas level-density models there is an implicit 
assumption of equal number of states with different par- 
ity at any given excitation energy. However, models that 
take into account microscopic effects show clear struc- 
ture in the parity ratio, see eg. Refs. [TH, H3, 5j|. In 
connection with astrophysical reaction rates the parity 
ratio may play an important role, and can be included in 
the Hauser Feshbach formalism, see eg Refs. [H. l45l] . 

From the single-particle point of view it is clear that 
the parity ratio should show structure. For ||Sr46 the sit- 
uation is illuminating. The nucleus is close to spherical in 
its ground-state (£2 = 0.05) [2l[. It has 2 proton holes in 
the pf shell and 6 neutrons in the <7g/ 2 shell. Since there 
are large gaps in the single-particle spectrum the energy 
to excite nucleons across the gaps, which cause changes 




Pairing Gap A [Me VI 

FIG. 3: (Color online) Proton pairing-gap distributions at 
different excitation energies 

£ cxc for 162 Dy shown in 20 keV 
pairing-gap bins. The proportion of paired states (A > 0) are 
shown in percent in the boxes. 
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FIG. 4: Average proton pairing gap (A) at different excitation 
energies -E oxc for 162 Dy. The mean values are given by the 
distributions in Fig. [3] 



in parity, is quite large. The parity ratio displays long- 
range oscillations, see upper left panel of Fig. O which are 
directly connected to the large single-particle gaps. The 
single-particle gaps effectively decrease with increasing 
deformation, and the oscillations are therefore expected 
to decrease with deformation. To test these ideas the 
parity ratio in 84 Sr is calculated at different deformations 
and shown in Fig. [5j Indeed, for larger deformations the 
oscillatory pattern vanishes and the parity ratio equili- 
brates to one. 

In Fig. [H the parity ratios for 56 Fe, 60 Ni and 68 Zn are 
shown and compared to two calculations of Ref. [l6[ ■ The 
solid blue lines show the CFY model, the red dashed 
lines show the statistical parity projection model and 
the black dots show calculations using the Shell-Model 
Monte-Carlo method. The parity ratio in the Fe and 
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FIG. 5: (Color online) Parity ratio versus excitation energy 
for 84 Sr calculated with different deformations. Large defor- 
mations imply a parity distribution close to 1 even at low 
excitation energies. For small deformations the parity non- 
equilibrium can survive to high excitation energies. 



FIG. 7: (Color online) Parity ratio versus excitation energy 
for Sr isotopes. The solid blue lines are obtained from the 
CFY model and the dashed red lines are obtained from the 
statistical parity projection of Ref. [45j . 




15 - Zn, e 2 =-0.15 
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FIG. 6: (Color online) Parity ratio versus excitation energy. 
The solid blue lines are calculated with the CFY model. The 
red dashed line and black dots are given b y a statistical model 
and a Monte-Carlo method, respectively [la ]. 



Ni isotopes is not equilibrated below 15 MeV in any of 
the calculations, while 68 Zn equilibrates at much lower 
energies (below 10 MeV). For 60 Ni and 68 Zn there is 
a clear oscillatory behavior prior to equilibration in the 
CFY model, and in the case of Ni there is a good agree- 
ment between the Monte-Carlo method and the CFY 
model, especially in the region of 8-12 MeV excitation 
energy. Note, however, that fluctuations seen in the 
present micro-canonical approach may be smeared out 
by the grand-canonical approach as used in Ref. [16j |. 

The parity ratio has been calculated for Sr-isotopes 
within a statistical method in Ref. [45j where the im- 
pact on astrophysical reaction rates was investigated and 
found to be small. The statistical method gives at most 
one oscillatory maximum before it equilibrates, as seen in 



Fig. [7] (dashed lines) . In the CFY model the parity ratio 
has substantially more structure (solid lines). Nuclei with 
small ground-state deformations show long-range oscil- 
lations which survive to high excitation energies before 
equilibration. The overall results are quite different from 
the model of Ref. [H- 

C. Angular momentum distribution 

In Fermi-gas models the distribution of angular mo- 
mentum is given by the Gaussian envelope in the spin 
cutoff model, see eg. Ref. Q, 

which is obtained by random coupling of uncorrelated 
spins of the nucleons. In Eq. (f2"2"|) the spin cutoff factor, 
<7, is defined by 

2 frigid U . . 

a= ~Va< (23) 

where is the rigid-body moment of inertia, U = 

E — S is the effective excitation energy shifted by the 
back-shift 5, and a is the level-density parameter. 

In Fig.[5]the angular-momentum distributions for 68 Zn 
and 162 Dy are shown for several excitation-energy re- 
gions. The black lines with dots show the CFY model 
results while the red solid lines show the Gaussian dis- 
tribution of Eq. (f!?2"|) with the spin cutoff factors given in 
Ref. y|, and the blue dashed lines show Gaussian distri- 
butions fitted to the combinatorial calculation. For low 
excitation energies there are clear deviations from the 
Gaussian profiles while for higher excitations the combi- 
natorial distribution tends to the Gaussian profile. How- 
ever, for 68 Zn the deviations from a Gaussian angular 
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FIG. 9: (Color online) The top panel shows the spin cutoff 
factor as a function of excitation energy for 162 Dy in the CFY 
model (black solid line) and the statistical model of Ref. [H 
(red dashed line). The bottom panel shows the ratio of the 
combinatorial and statistical spin cutoff factors of the top 
panel. 



the CFY model. 



IV. COMPARISON WITH EXPERIMENTAL 
DATA 

A. Neutron separation-energy level spacings 



FIG. 8: (Color online) Angular-momentum distribution for 
68 Zn and 162 Dy in 4 different excitation-energy regions. The 
black lines with dots show results from the CFY model. The 
red solid lines and blue dashed lines show the Gaussians given 
by Eq. (|22j> using the spin cutoff factors from Ref. [4| and a 
direct fit to CFY, respectively. 



momentum distribution remain up to about the neutron 
separation energy (5„=6.48 MeV). 

For low excitation energies in 68 Zn there is an odd-even 
spin staggering which is not explained by the spin cutoff 
model. This effect has also been observed in Fe-isotopes 
in the Shell Model Monte-Carlo Method 

Fig. [H] shows the spin cutoff factor deduced from the 
CFY model (black solid line) and from the statistical 
model [H (red dashed line), as a function of excitation 
energy for 162 Dy. For this particular nucleus the spin 
cutoff factor for excitation energies > 3 MeV is similar in 
shape but ~lfj % larger than in the statistical model 
The fact that the statistical model is smaller than the 
CFY model, despite the presence of pairing in the CFY 
model, is accidental for this nucleus. Since both the rigid 
body moment of inertia and the level density parameter 
are fitted to experiments using only three parameter, the 
detailed description for a particular nucleus might give 
this result when compared to a very different model as 



The s-wave neutron resonance spacings constitute the 
most comprehensive experimental database for compar- 
ison with NLD calculations [4a|. This database serves 
as a benchmark for all large-scale level-density models 

The s-wave neutron resonance spacing Dq at the 
neutron separation energy S n of the compound nu- 
cleus (Z,N) is obtained from calculated level densities 
p(E,I,n) as, 

1 

Do ~ 

P (S n , h + 1/2, 7T ) + p (S n , I - 1/2, 7T ) for I > 0, 

p(S n ,l/2,7r ) for I -0, 

(24) 

where Iq is the ground-state spin and ttq is the ground- 
state parity of the target nucleus (Z, N — 1). In Fig. [TU] 
we study the ratio of calculated and experimental level 
spacings at the neutron separation energy, -DTh/£>Ex P , 
for all nuclei in the database. 

The quality of a level-density model can be estimated 
by the rms-factor 



/rms — exp 



N e 



, ln 2 -^ 

i=l "^Exp 



n 1/2 



(25) 
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FIG. 10: Ratio between theoretical and experimental level 
spacings at neutron separation energy versus mass number 
A. The rms-factor is / rm s = 4.2. The experimental data are 
taken from Ref. [4611. 



and the mean factor 



exp 



N, 



N r ^ 



D 



Th 



D 



Exp 



(26) 



where D^ h and -Dfj xp are the theoretical and experimen- 
tal level spacings and N e is the number of nuclei in the 
database. The CFY model gives / rms = 4.2 and m = 1.1, 
see Fig. [10] of the error ratio versus A. These results are 
compared to other statistical and combinatorial models 
in Table HT1 and discussed in Sec. [Vl 

There seems to be a clear residual shell structure in the 
level spacings, as seen in Fig.[TUl especially in the doubly 
magic 208 Pb region. A similar effect can be observed 
for the Gogny model in Ref. [15| . and in results of the 
Skyrme model with the BSk9 interaction, as shown in 
Ref. [14(. In addition, there seems to be a drift with 
mass number where especially nuclei in the rare-earth 
and actinide regions over-estimate the neutron separation 
level spacings. A similar drift can be seen in the Gogny 
model of Ref. [l5| while it seems not to be present in the 
Skyrme-HFB model in Ref. [3]. 

The rare-earth and actinide regions, where too large 
level spacings are calculated, mainly correspond to de- 
formed nuclei. On the other hand, nuclei corresponding 
to regions with particularly low ratio -Dxh/^Exp seen in 
Fig. [10] (around ^4=32, 132 and 208) correspond to spher- 
ical nuclei. The overestimation of D for deformed nuclei 
(too low calculated level density), and underestimation 
of D for spherical nuclei (too high calculated level den- 
sity) seems to be approximatively systematic, as seen in 
Fig. 1111 where the ratio between theoretical and experi- 
mental level spacings is shown versus the absolute value 
of the deformation. 

However, as seen in the lower panel of Fig. [TTJ there 
is no clear correlation with the microscopic energy E m { c 
of Ref. [2l[ as might be expected from the residual shell 
structure seen in Fig. [TOl 



FIG. 11: (Color online) The top panel shows the ratio between 
theoretical and experimental level spacings at neutron sepa- 
ration energy versus the absolute value \e2\ of the calculated 
quadrupole deformation. The bottom panel shows this ratio 
as a function of the microscopic energy. The solid red line 
shows an exponential fit to illustrate the correlation. The mi- 
croscopic energies and the quadrupole deformations are taken 
from Ref. [2lJ. The experimental data are from Ref. (4f|. 



B. Detailed level-density functions 

The Oslo method is a commonly used experimental 
method for extracting detailed level-density functions for 
large ranges of excitation energy [19( . It provides a valu- 
able test for NLD models. However, the approach is 
model dependent since the level density is extracted by 
use of a back-shifted Fermi Gas approximation and a 
spin-cutoff factor, see also discussion in Ref. [38j]. 

In Figs. [T2] and [T3] level densities obtained from the 
CFY model are compared to data for a number of nuclei 
where experimental data are available [48|, |4tJ [5(| HH [H, 
El. 

In Fig. [T2]we show level densities of the rare-earth nu- 
clei 148 ' 149 Sm, 161,162^ 166,167^ and 170,171,172^ as 

functions of excitation energy. The data are in general 
well reproduced by the CFY model, with an error of less 
than a factor of 2. The good agreement between the cal- 
culated and experimental slopes indicates that the single- 
particle structure and the moments of inertia in the rota- 
tional bands are sound. However, an observable trend in 
this mass region is that the level density for even-even nu- 
clei is slightly over-estimated and for odd nuclei the level 
density is slightly under-estimated. The effective back- 
shift is to a large extent controlled by the ground-state 
pairing gap. By fine-tuning the pairing gaps it is possi- 
ble to get an almost perfect agreement with experiment. 
However, in this paper no efforts are made to adjust pa- 
rameters of the model to fit measured level-density data. 
Indeed, the ground-state pairing gaps are given by the 
mass model, see Sec. IIIBI Neither are any renormaliza- 
tions of calculated level densities to fit data performed, 
as is sometimes done in other models, see Refs. [lH l38j. 
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FIG. 12: (Color online) Level densities p as functions of ex- 
citation energy for Sm, Dy, Er and Yb isotopes. The black 
solid lines show the CFY model and the red dots show the 
experimental data [Hoi. l5ll. IH^. . 
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FIG. 13: (Color online) Level densities as functions of excita- 
tion energy for 50,51 V and 97,98 Mo. The black solid lines show 
the CFY model and the red dots show the experimental data 

mm. 



ing, see Fig. [TOJ which have an overall / rms = 4.2. Also, 
the slopes and the detailed structures are not as well de- 
scribed in these nuclei as in the rare-earths. Especially 
50 V exhibits an oscillatory pattern in the CFY model, 
which is an effect of the small ground-state deformation 
£2 = 0.05 2l|. A corresponding oscillatory pattern is not 
as clearly present in the experimental data. In 51 V the 



experimental data shows oscillatory behavior while the 
CFY model is more smooth. The ground-state deforma- 
tion is slightly larger with £2 = 0.083 [2l[, which in the 
CFY model means that since the moment of inertia is 
larger the rotational states have a larger influence and 
smooths the level density. 



Experimental level-density functions are available also 
for lighter mass regions, and in Fig. [T3] data for V and 
Mo isotopes are shown. The overall agreement between 
the model and these experimental data is somewhat in- 
ferior to what was obtained in the rare-earth region. For 
the Mo isotopes in Fig. Q2] the CFY model is roughly a 
factor of 3 too large for high excitation energies while 
for the V isotopes the over-estimation is roughly a fac- 
tor of 4. However, these errors are consistent with the 
overall results from the neutron separation level spac- 



C. Low- lying discrete levels 

The amount of experimental information on low-lying 
discrete energy levels is indeed substantial. In addition 
to the energy parity and angular momentum are often 
known. In Ref. [13] a global comparison is made between 
measured and FRLDM calculated ground-state spins and 
parities for odd-mass nuclei, and the agreement was 
found to be very good. The collection of low-energy data 
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provides important information for level-density models, 
since the accuracy of the model can be studied at the 
very lowest energies. Cumulative distributions of known 
levels can thus be constructed, and compared to calcula- 
tions. In Ref. [lj] a selection of 15 different nuclei was 
made to represent all different kinds of nuclear aspects, 
such as light-heavy, spherical-transitional-deformed, as 
well as odd-odd, even-even and odd-even nuclei (also used 
in Ref. [HI). We therefore compare our model low- lying 
level densities to experimental data for the identical set 
of 15 nuclei. 

The cumulative level density for the 15 nuclei are 
shown in Fig. [T5J The CFY model (solid blue line) is 
compared to experimental data (black dash-dotted line) 
and to the HFB model (red dashed line) of Ref. [Hj]. For 
a good agreement with data, the calculated line should 
follow as precise as possible the experimental curve at 
the lowest energies, and should then smoothly deviate, 
always larger than data. 

In general, the theoretical models seem to give very 
comparable results for deformed nuclei while differences 
are larger for transitional and spherical nuclei. For the 
light spherical nuclei, like 42 K and 56 Fe there are pro- 
nounced single-particle gap effects which cause jumps in 
the level density. These effects are smoothly disappear- 
ing at higher energies. A quite drastic deviation between 
our calculation and data is seen for 208 Pb. The calcu- 
lated (cumulative) level density is systematically larger 
starting already at about 3 MeV excitation energy (in 
the figure seen at 0.5 MeV; the energy scale is shifted by 
2.5 MeV for this nucleus). This is also seen in the to- 
tal level density at the neutron separation energy, and in 
Table U the ratio between theoretical and measured level 
spacings at the neutron separation energy is listed for the 
9 nuclei. For 208 Pb this ratio is 0.02, i.e. the calculated 
level density is about 50 times larger than experiment! 

When com par ing to the HFB model calculation of the 
level density |38j one should note that interpolations be- 
tween different deformations are included in this model, 
while in our model the ground-state deformation is kept 
fixed for all excited states. In particular, this somewhat 
phenomenological way to treat deformation changes has 
strong effects on the results for 127 Te. 
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FIG. 14: (Color online) Accumulated level density as a func- 
tion of excitation energy. The blue solid and red dashed lines 
show the present CFY model and the HFB model [H- The 
black dash-dotted line shows experimental data collected in 
Ref. 



D. Parity ratio 

The parity ratio has been measured experimentally by 
Kalmykov et.al. [2(| for the two spherical nuclei 58 Ni 
and 90 Zr. For these nuclei the 1 — 2 angular- momentum 
component of the level density is measured and separated 
into parity components. In the top panel of Fig 
the level-density component p(E, I = 2, tt = ±1) for 
is compared to CFY model calculations. The level den 
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TABLE I: Table of neutron separation energy spacings com- 
pared to experimental data for the nuclei showed in Fig. 1141 



sities from the CFY model are in good agreement with 
experimental data, and in the lower panel of Fig. 15(a 
the parity ratio is shown. Predictions for the total level 
density as well as the 1 = 2 component are shown as black 
solid and red dashed lines, respectively. The difference 
between the parity ratio for the total level density and the 
1 = 2 component decreases with excitation energy, be- 
cause the spin cutoff model becomes more realistic when 
the excitation energy increases, see Sec. IIII CI For the 
parity ratio the CFY model result is within the experi- 
mental error-bars for all excitation energies and seems to 
show a similar pattern as the experiments: a high par- 
ity ratio (excess of negative parity states) at 8 MeV and 
a low ratio (excess of positive parity states) at 11 MeV 
excitation energy. The blue dot-dashed line shows the 
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FIG. 15: (Color online) The top panels show CFY calculated 
components of the 2 + (black solid) and 2~ (red dashed) level 
densities for 90 Zr (top part of figure) and 58 Ni (lower part 
of figure), as functions of excitation energy. Experimental 
data from Ref. ^20] are shown as black dots and red diamonds 
with error-bars for 2 + and 2~, respectively. The lower panels 
show for each nucleus the parity ratio versus excitation energy. 
The solid black and dashed red lines show the CFY model 
results for the total level density and the 1=2 component, 
respectively. The blue dot-dashed lines show the Skyrme- 
HFB model of Ref. [JJ|]. Experimental data from Ref. [20] are 
shown with error-bars. 



Skyrme-HFB model of Ref. 14]. It is also within the 
experimental error-bars for all excitation energies except 
at 9 MeV where the HFB model gives a large positive 
ratio (~ 3) while the experiments and the CFY model 
are close to unity. 

The 1 = 2 component of the level density in 58 Ni is 
shown in the top panel of Fig. |15(b)| The agreement 
between experimental data and the CFY model is some- 
what inferior to what we obtained for 90 Zr. For both 
parities the level density seems to increase faster in the 
CFY model than what is seen in the experimental data. 
At 14 MeV excitation energy the model overestimates 
the level density by roughly a factor of 6. The parity 
ratio is shown in the lower panel of Fig. 15(b) The ex- 



T3 

"S io s 



90 Nb f=l + 



5 6 7 8 9 10 

Excitation Energy [MeV] 

FIG. 16: (Color online) Level density of the 1 + component of 
the level density as a function of excitation energy for 90 Nb. 
The black solid and red dashed lines show the CFY model 
and Skyrme-HFB model of Ref. respectively. Data from 
Ref. [20|] are shown as black dots with error-bars. 



perimental parity ratio is close to unity at 8 MeV and 
then decreases with increasing excitation energy. The 
CFY model shows a different pattern. It gives a very 
low parity ratio at low excitation energies with an al- 
most monotonic increase with increasing excitation en- 
ergy. The ratio only becomes close to unity at 20 MeV. 
The Skyrme-HFB model of Ref. [l4| shows a similar pat- 
tern as the CFY model but with an even lower ratio for 
excitation energies below 20 MeV. 

The 1 + component of the level density of 90 Nb has 
also been measured in Ref. 20]. Fig. [TBI shows the ex- 
perimental data compared to the CFY and Skyrme-HFB 
models. In the experimental data there is a clear os- 
cillating structure. The CFY model (black solid line) 
over-estimates the level density and the oscillating struc- 
ture is much less pronounced. In the Skyrme-HFB model 
(red dashed line) there are long-range oscillations simi- 
lar to what is seen in the experimental data, but the 
energy separation between consecutive minima is larger. 
The level density is under-estimated in the Skyrme-HFB 
model, whereas the CFY model over-estimates the level 
density by a similar factor. 



V. COMPARISON WITH OTHER MODELS 

The CFY model is here compared to other statistical 
and combinatorial NLD models that provide neutron res- 
onance level spacings. In general, the statistical models 
listed in Table HII have an rms deviation just below 2. This 
low rms deviation is obtained because several parameters 
in the level-density formulas are directly adjusted to the 
neutron separation-energy level spacings and low-lying 
discrete energy levels. 

The back-shifted Fermi model of Ref. (4?| is based on 
a simple Fermi-gas formula whereas the Constant Tern- 
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Statistical Models 


/rms 


Ref. 


Back-shifted Fermi Model 


1.71 


[4Jj 


Const. Temp. Model 


1.77 


[47] 


Back-shifted Fermi + Const. Temp. 


Model 1.7 


1 


Generalized Superfluid Model 


1.94 


[47] 




Combinatorial Models 


frms 


Ref. 


Skyrme-HFB 


2.35 


[47] 


Gogny-HFB 


4.55 


[15] 


CFY 


4.2 


Present 



TABLE II: Table of rms-factors / rms for statistical and com- 
binatorial models for neutron separation energy level spacings 
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FIG. 17: Vibrational enhancements at the neutron separation 
energy in the attenuated phonon method (top panel) and the 
QTDA method (bottom panel) versus of mass number A. 



perature model is based on the approach of Gilbert and 
Cameron [10| . Both models contain explicit enhance- 
ment factors for rotations and vibrations. The model 
by Rauscher, Thielemann and Kratz 0] is also based on 
the approach of Gilbert and Cameron. In contrast to 
the first two models this model has no explicit enhance- 
ment factors for rotations and vibrations. It accounts for 
shell effects and thermal damping in terms of an effective 
level-density parameter. The model uses the microscopic 
energy corrections of Ref. [2l[ together with three free 
parameters in the fitting procedure. The Generalized Su- 
perfluid Model is similar to the models mentioned above. 
In addition it takes into account how pairing evolves with 
increasing excitation energy. It also incor por ates explicit 
rotational and vibrational enhancements [471 ] . 

Predictions by statistical models in regions outside 
the fitting region probably give much less accurate re- 
sults than in the adjustment region. On the other hand, 
since the combinatorial models are all based on calcu- 
lated single-particle spectra they could in principle have 
better predictive power, in particular in regions where 



the single-particle model is sound. In contrast, they are 
not equally flexible in terms of parameter fits to the level 
spacings database. Therefore, the rms deviation factor 
of combinatorial models with respect to known data is 
larger than for the statistical models. The / rms = 4.2 
for the CFY model is about a factor of 2 larger than in 
statistical models. However, the latter result is obtained 
with no parameters specifically fitted to the level-density 
data. The mean deviation factor m = 1.1 indicates that 
the level density is on average reasonable. However, there 
seems to be a drift with mass number, as seen in Fig. [TQl 
present in the CFY model which is not seen in the HFB 
models of Refs. [II, ll}. 

In the lower part of Table ID the CFY model is com- 
pared to two other large-scale combinatorial NLD mod- 
els based on the HFB method. The model of Ref. [l5[ 
is based on the Gogny D1S interaction for the mean- 
field and incorporates combinatorial rotations (similar to 
Sec. Ill CI but no pairing dependence of the moment of 
inertia), and the attenuated phonon method is applied 
to account for vibrational states. Pairing is included 
explicitly for the ground-state, and excited states are 
back-shifted by an energy-dependent gap procedure. The 
model gives / rms = 4.55 for the subset of even-even axi- 
ally deformed nuclei. The error is slightly larger than the 
CFY model. 

For the Skyrme-HFB NLD model the deviation factor 
is /mis = 2.35 [l?]], which is ~ 35% larger than the sta- 
tistical models and 44% smaller that in the CFY model. 
This model is based on a Skyrme-HFB mean-field to- 
gether with combinatorial rotations and the attenuated 
phonon method for modeling vibrational states. In ad- 
dition a phenomenological deformation change from de- 
formed to spherical shape at some specified excitation 
energy is incorporated. 

The attenuated phonon method is a phenomen olog ical 
way to model nuclear vibrations, see eg. Refs. [HjEj]- It 
assumes that the quadrupole and octupole phonon states 
in nuclei can be modeled by a gas of non-interacting 
bosons. The vibrational excitation energies are taken 
from systematics of the lowest non-rotational 2 + and 3~ 
states. The model is formulated as a multiplicative factor 
^vibr = exp [SS — SU/T) which enhances the level den- 
sity. T is the nuclear temperature and the SS and SU are 
the entropy and internal energy change induced by the 
bosons. The occupation probabilities of the bosons are 
described by damped Bose statistics, where the damping 
sets in at considerably higher energies than the neutron 
separation energy. At excitation energies close to or be- 
low the neutron separation energy the damping has neg- 
ligible effect which implies that the phonons are allowed 
to be repeated several times. 

The vibrational enhancement at the neutron separa- 
tion energy is quite small in the QTDA method as com- 
pared to the attenuated phonon method, see Fig. [TTj 
The QTDA gives an enhancement of the order of a few 
percent compared to up to a factor 7 in the attenuated 
phonon method. The largest vibrational enhancements 



14 



are obtained for the A=75 region and for Cd isotopes in 
the A=115 region, which is consistent with the strong 
(quadrupole) vibrational character of nuclei in these re- 
gions. 

The QTDA phonons are microscopically built from 
quasi-particle excitations E^' 1 with energies not much 
different than the QTDA phonon energies (hoj)j- The 
way to fully account for double-counting of phonon states 
(Eq. I19|) . thus implies a small vibrational enhancement 
since the existence of a few low-lying collective phonons 
will not impact the level density if all other phonons 
are very non-collective. The non-collective character of 
most of the phonons also implies that they cannot be re- 
peated. Phonons in the attenuated phonon model are not 
described microscopically, and double-counting of states 
may clearly appear. In addition, it assumes that each 
phonon can be repeated several times. 

VI. SUMMARY 

A combinatorial model for the nuclear level density 
is presented. The model is based on the folded- Yukawa 
single-particle model with ground-state deformations and 
parameters from Ref. [2lj]. The model is used to calcu- 
late the neutron resonance level spacings, yielding an rms 
deviation of / rms = 4.2. It also compares favorably with 
experimental level density data versus excitation energy 
for several nuclei in the rare-earth region as measured 
by the Oslo method, as well as with the cumulative level 
density extracted from low-energy spectra. 

The role of collective enhancements has been investi- 
gated in detail. Pairing is incorporated for each individ- 
ual many-body configuration, and the distribution of the 
pairing gaps is investigated. No sharp pairing phase tran- 
sition is observed. Instead, even at the highest excitation 
energy considered, a non-negligible fraction of the states 
have a considerable pairing gap. 



Rotational states are included combinatorially by a 
simple rotor model with a moment of inertia dependent 
on the deformation and pairing gap. Vibrational states 
are included using a Quasi-particle- Tamm-Dancoff Ap- 
proximation. It is found that the vibrational enhance- 
ment in the QTDA model is very small, on the order of 
a few percent at the neutron separation energy. 

The parity distribution in the CFY model shows large 
oscillatory patterns for nuclei which have large gaps in 
the single-particle spectrum, separating shells with dif- 
ferent parities. This is often the case for nuclei with 
small deformations. The patterns are quite different from 
the smooth pattern of the statistical parity distribution 
model of Ref. [l6[ . The CFY model is compared to other 
models and to experimental data when available. 

The angular-momentum distribution in the CFY 
model is compared to the spin cutoff model in Sec. IIII CI 
It is found that the Gaussian envelope of the spin cutoff 
model is in good agreement with the CFY model for high 
excitation energies. 
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